o 
o 

(N 



(N 



O 

o 



Isotropic cosmological singularities in 
spatially-homogeneous models with a cosmological 

constant 



Paul Tod* 

Mathematical Institute and St John's College 
Oxford 

< 

q\ ! February 1, 2008 



cr, Abstract 

j^jrj' We prove well-posedness of the initial value problem for the Einstein 

equations for spatially-homogeneous cosmologies with data at an isotropic 
cosmological singularity, for which the matter content is either a cosmo- 
logical constant with collisionless particles of a single mass (possibly zero) 
' or a cosmological constant with a perfect fluid having the radiation equa- 

tion of state. In both cases, with a positive cosmological constant, these 
solutions, except possibly for Bianchi-type-IX, will expand forever, and be 
geodesically-complete into the future. 



PACS: 04.20EX, 98.80Jk 



1 Introduction 



• i 

An isotropic cosmological singularity is one which may be removed by confor- 
mally rescaling the metric. More accurately, the physical space-time (M,g) is 
conformally related to an unphysical space-time (M, g) by g = ft 2 g where f2 is 
a function positive on M and vanishing on a smooth space-like hypersurface £ 
in M, which provides a boundary to M at which g is regular but g is singular. 

In several cases of interest, it is possible to formulate a well-posed initial 
value problem in M for the Einstein equations in M by giving data at £. 
This was done for a perfect-fluid source with a linear equation of state in [1], 
for massless Einstein- Vlasov with spatial homogeneity in [2] and for massless 
Einstein- Vlasov without the assumption of spatial homogeneity in [3]. Steps 
towards the same end for a range of other matter models were taken in [4] and 
our purpose here is to carry through the proof of well-posedness in two cases 
with a cosmological constant A, namely massive or massless Einstein- Vlasov, 
and perfect fluid with the radiation equation of state. 

One motivation for this calculation is the recent appearance of long-time 
existence results for massive Einstein- Vlasov and a positive A [5], [6]. In [5], 
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it was shown that spatially-homogeneous cosmologies of Bianchi-types other 
than type- IX, once expanding, would expand forever and would be geodesically- 
complete into the future. (Type- IX may do this but could recollapse instead.) 
In [6], the corresponding result was proved for plane or hyperbolic symmetry. 
Putting what we find here together with Lee's result in [5], we may construct 
solutions of the Einstein equations with this matter model or with massless 
Einstein- Vlasov which expand away forever from an isotropic cosmological sin- 
gularity to a future infinity or Scri (in the sense that the metric can be rescaled 
to add a future boundary at infinity; strictly speaking, to call this Scri, we 
should check for appropriate differentiability of the rescaled metric but, for the 
Einstein- Vlasov cases, we don't do that here). Given that result, it is natural to 
look for perfect fluid cosmologies with positive A which also expand away forever 
from an isotropic cosmological singularity towards a Scri (this time, appropri- 
ately differentiable) , and that is the other concern in this article. Both classes 
of models then provide fairly explicit examples of the cosmologies considered 
recently by Penrose [7] : the conformal metric can be extended both through the 
initial singularity (or bang) and through the future infinity (or Scri). However, 
note that Penrose wanted more than this, namely that the conformal geometry 
at bang and Scri were essentially equivalent, and these models don't seem to 
have that property, though this isn't completely clear. 

For the Einstein- Vlasov equations, our technique is simply to re-examine 
the proof of [2], with massive Vlasov in place of massless, and including a 
cosmological constant. As one would expect heuristically, and as was suggested 
without complete proof in [3J, these modifications are unimportant near the 
singularity. Our result is Theorem 12.11 of section 12.31 that there is a well-posed 
initial value problem with data at the singularity, and the data is just the initial 
distribution function, subject to the vanishing dipole condition as in [2]. 

As Corollary 12.11 we prove local-in-time existence for massless Einstein- 
Vlasov with A near an isotropic singularity, and then for completeness prove 
long-time existence for this matter model for positive A (Theorem I2.2|) . 

For the perfect fluid with radiation equation of state and A, we follow the 
calculation of [lj to obtain local existence in time for a variety of spatially- 
homogeneous cosmologies, and then directly prove long-time existence. 

The plan of the paper is as follows. Section 2 is concerned with the Einstein- 
Vlasov results and Section 3 with the radiation results. 

We begin in sub-section 2.1 by collecting together some results on conformal 
rescaling for the massive Einstein- Vlasov equations with cosmological constant. 
In subsection 2.2, we specialise to the spatially- homogeneous case and write 
the Einstein equations as a first-order system of ordinary differential equations, 
coupled to the Vlasov equation which is still a partial differential equation. We 
discuss the splitting of the Einstein equations into constraints and evolution, 
and the result of imposing the constraints at the singularity. In subsection 2.3, 
we state the existence theorem, and show how it can be proved by the same 
techniques as in [2] for the massless Einstein- Vlasov system without A. The 
proof can then be seen to carry through for massless Einstein- Vlasov with A. In 
subsection 2.4, closely following [5], we prove long-time existence for massless 
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Einstein- Vlasov with positive A, so that for both the massive and massless case 
there exist solutions of the kind imagined by Penrose, with the conformal metric 
extendible through both the bang and Scri. However, it seems that even in the 
conformal geometry, the bang and Scri can be distinguished. 

In Section 3, for the radiation results, we begin in subsection 3.1 with the 
Friedman-Robertson- Walker (FRW) metric. In terms of conformal time, the 
Einstein equations reduce to the Friedman equation for the scale factor. The 
scale factor grows from a simple zero at the bang to a simple pole at Scri, after 
a finite conformal time. For positive A, the Friedman equation has a remark- 
able symmetry in that the inverse of the scale factor satisfies the same equation 
as the scale factor, but with the density parameter and cosmological constant 
interchanged. This transformation also interchanges bang and Scri but, as we 
shall see, its existence depends on the absence of conformal curvature. In sub- 
section 3.2, we consider the Bianchi type-Ill and Kantowski-Sachs metrics. We 
show that there are solutions expanding from an isotropic singularity at the 
bang. For positive A, the Kantowski-Sachs metrics may recollapse but the 
Bianchi-III metrics cannot and must expand to Scri, again after a finite con- 
formal time. Finally, and most generally, in subsection 3.3 we consider metrics 
of Class A in the Bianchi classification. All types admit solutions expanding 
from an isotropic singularity at the bang. For positive A, all but type- IX, which 
may recollapse, expand to Scri, again after a finite conformal-time. The sym- 
metry seen in FRW does not persist in the other examples: both the bang and 
Scri are umbilic in the conformal metric but the conformal metric is even in 
conformal-time about the bang but apparently not about Scri. 

2 The Einstein-Vlasov equations 

2.1 Conformally rescaling the massive Einstein-Vlasov equa- 
tions 

The physical space-time will be (M,g), while the rescaled, unphysical space- 
time will be (M, g) with 

g ab = ft 2 g a b- (1) 

The physical space-time has an isotropic cosmological singularity if vanishes 
at a smooth, space-like hypersurface £ in the unphysical space-time M. 

In kinetic theory, the matter content of space-time is taken to consist of 
a collection of particles which move on geodesies between collisions. States 
having a given mass m > define thejseven-dimensional submanifold P m of 
the one-particle phase-space P = T*M. On P m we take (x a ,pi) as local 
coordinates, po being determined by the equation g ab p a Pb = m 2 , an d the 
requirement that p a be future directed (conventionally for us, a, b, c, .. = — 3 
while i,j,k,.. = 1 — 3). It is important for conformal rescaling to retain the 
notational convention p a = g ab pb- 

The free-fall trajectories define on P a congruence of curves, along which 
H = \g ah p a Pb is constant. The geodesic spray L~ g is the vector field tangent 
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to these curves defined, in local coordinates, by 



d 1 dg ab 
" ' Pa dxb ~ 2 PaPb ^ d Pc 



^9=~9 ab Pai- h -^PaP b ^-r^-. (2) 



The cotangent space to M at a point x is a flat Lorentzian manifold, and 
on the submanifold P m (x) there exists an invariant volume measure Cb m given 
by 

1 „ 

Cu m = ——==dpi. (3) 
PW-9 

In this article we shall be concerned with particles of a single mass, so that we 
shall henceforth omit the subscript m on Co. The distribution of particles and 
momenta is described by a non-negative scalar function / = f(x a ,pi) on P m , 
and the condition that / represents a collisionless gas is simply 

W) = 0, (4) 

which is conveniently called the Vlasov equation. Note that / satisfies the 
Vlasov equation if and only if it is constant along geodesies of g. 

The stress-energy-momentum tensor due to these particles of mass m is 
given by 



Tab(x) = / fPaPb (5) 

JPm{x) 

and if the Vlasov equation is satisfied then 

V a T a6 = 0, (6) 

where V is the metric covariant derivative for g. 

The coupled Einstein- Vlasov equations, for the metric g a b and the particle 
distribution function / are therefore 



Gab = 8vr / fp a p b Co 

J Pm 

W) = 

taking G = c = 1. 

We shall investigate the transformation properties of the Einstein- Vlasov 
system under conformal rescaling. Under rescaling we choose, as our notation 
already implies, 

Pa = Pa, 

so that the canonical one-form is unchanged: 

6 = p a dx a = p a dx a . 

This then necessitates 

g ab p aPb = n 2 g ab PaPb = tfm 2 , (7) 
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and at the approach to the initial surface S in M (which is the singularity in 
M), where O vanishes, the particles are effectively massless. Next 

to / 5 

OPa 



c § f = n-\£ g f + m 2 n 2 r a - 



where, as usual, T a = f2 d a Q, and the distribution function / does not change 
under rescaling. Thus, in M, the Vlasov equation from M becomes the equation 



C g f + m 2 fi 2 T a |^ = 0. (8) 

dp a 



Again, the mass term vanishes at S. For the volume- form on the mass-shell we 
obtain 



1 

p u V~-g 

and we can define a rescaled energy-momentum tensor as 



T ab = tt 2 T ab = / fpaPb to 



p 

This will not be divergence- free, since in fact 

~g ab V a T bc = Q-\g ab V a T bc - g ah T ab T c ), 

where V is the metric covariant derivative for g. Here the left-hand-side van- 
ishes, which gives an equation for the divergence of T ab from the right-hand-side. 
Note the following equation for the trace: 

T = g ab T ab = n 2 m 2 j fuj. 

As one of the conditions for an isotropic cosmological singularity, we assume 
that / extends to a smooth function on T*£. Therefore the trace of the rescaled 
energy-momentum tensor vanishes at S, which is a reflection of the earlier 
observation that the mass of the particles is negligible there. 

With the conformal transformation of the Einstein tensor taken e.g. from 
PP, we can write down the conformal EV equations for g ab and / as 

n 2 G a b = 2W a v b n- 4v a nv b n- g ab {2Qnn- v c nv c n) 

[ fpaPb^lr + An* 9ab (9) 
V^g J p u 

together with (jSJ). 

2.2 The spatially-homogeneous case 

We briefly recall the formalism for spatially-homogeneous cosmologies (see e.g. [8]) 
These have a three-dimensional isometry group G transitive on space-like hy- 
persurfaces. Suppose G has structure constants C l - k , for i,j,k = 1,2,3, then 
one introduces a basis of left-invariant one-forms a 1 satisfying 



^ = lcV^Aa fe . (10) 
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The metric is given in this basis by the expression: 

g ab dx a dx b = dt 2 - a„(i)oV (11) 

where t is proper-time on the congruence orthogonal to the surfaces of homo- 
geneity (which is necessarily geodesic) and dij is the physical metric on the 
surfaces of homogeneity. 

The Einstein- Vlasov equations reduce to a system of ordinary differential 
equations for d^, coupled with the Vlasov equation, which is a partial differen- 
tial equation. We seek a well-posed formulation of these in the rescaled space- 
time. Following the massless case from [2], we introduce a new time-coordinate 
r by t = t 2 and choose the conformal factor O = r. (The justification for 
making these choices is that they lead to a well-posed initial value problem.) 
The rescaled metric is then 

g ab dx a dx b = fl~ 2 g ab dx a dx b 

= ^-fljjMffV', (12) 

where a« = t~ 1 dij, which we shall assume is regular at the initial surface £ 
where r = 0. We shall systematically use b %3 for the matrix inverse to aij. 
We write the momentum in the basis of invariant one-forms as 

p = 2p dT + Pi(T % , 

and p a = ipQiPi). The normalisation ([7]) implies 

(pof = (P°) 2 = b ij PlPj +m 2 T 2 . (13) 

(Since p° = po with these conventions, we shall consistently use the form po 
to avoid later confusion with initial values.) For the Vlasov equation, we first 
obtain the geodesic equation, which is: 

d Pi a b 

dX = 7 * PaP ' 

where A is proper time for g, and the Ricci rotation coefficients j a bc are given 
by 

7% = = J a b0 

and 

i }_y ,m (C n n l C n n _|_ C n n \ 

Ijk ~ 2 ^ jm a nk t o km u nj t <-> jk a nm)- 

The geodesic equation can therefore be written 

^ = -C%V™p kPm . 

The apparent sign-change here is because the spatial metric, with our conven- 
tions, is —atj. However, we shall henceforth observe the convention that indices 
. . . will be raised with W and lowered with a«. 

The Vlasov equation simplifies in the spatially-homogeneous case because 
0, is a function only of time, so that T has only a zero-component and the 
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extra term drops out of ([8]). We assume that the distribution function will be 
homogeneous and distinguish the distribution function as a function of proper 
time t from the distribution function as a function of conformal time r by writ- 
ing f(t,pi) = f(r,Pi). For the homogeneous distribution function the Vlasov 
equation reduces to the statement that / is constant along geodesies of g a {,: 



# = 1 d£_d£ 
dX 2 dr dpi 



/>n— - ^-C k ijPk pib lj = 0, 



or 

^ = 2(b mn PmPn + m *T 2 y 1/2 C%PkPlb lj ^- i . (14) 

Now we follow |2j to find a first-order form of the equations, introducing a 
tensor fcj,- proportional to the second-fundamental form: 



d_ 



hi (15) 



jL b a = -b im V n k nm (16) 
ar 

— k{j — &Rij ~\~ [%Z%j 2kij b kmn&ij^) ~\~ kimkjnb 

-\^ mn k mn )k l3 + 8Ar 2 aij + ^m 2 a i3 [ ^ (17) 
2 v a J Po 

where Rij is the spatial Ricci tensor, so that 

SRij = -4:C k ck (C r tj a ir + C r ti aj r )b ct 
-AC c kl (C k c ,+C m tl a cm b kt ) 



and 



^ = K^/ /ftp 4f _ay )- (18) 

We have introduced the new variables Z%j, following [2], so that the pole in r in 
(|17p is no worse than first-order, but we therefore need an evolution equation 
for Zij. From the definition (118p . using (|14p . we find 



^ ry 1 t rr 7 647F 1 In f 9 f ^P 



' ° ;T f (( Po ) 2 b mn k mn + b m % nr k qrPmPn + 2m 2 r)fp iPj - d " P 



bo) 3J ' 

which again has no worse than a first-order pole. For later use, we write (|19p 
in the form 

— ~( ^ij kij I%j e/jj k mn ) Lij, (20) 
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where 



647r i. f df (ftp . . 

* = -Ta C - h Jdt PiPjPm ^ (21) 

J™ = ^| J f(( Po )H™ + b m «b nr p q p r )p lPj -^, (22) 

32vrm 2 r /" d 3 p 



From (JTSD and (jT6|) . 

-^-(a^ m - <5D = -k rn b rm ( aij V n - Sf), 

CLT 

so that, by a Gronwall estimate, if a^V™ = 5™ at r = 0, then this holds for all 
r. 

For the Hamiltonian constraint we calculate: 

C := ^ / / P0 d 3 p - | - r 2 (i? - -U«fcy + ±k 2 ) - \kr + 2Ar 4 = 0, (24) 
J 2 lo lb 2 

while for the momentum constraint we calculate: 

d:=^J f Pi d 3 p - r 2 (b mn k nj C j mi + b mn k ni C j mj ) = 0. (25) 

We calculate for the evolution of the constraints: 

d T (T 2 aC) = Cib^C m jm , 
d T (^d) = 0. 

so that the constraints are satisfied at all times if satisfied initially. 

We use a superscript to denote the value of a quantity initially. To satisfy 
the constraints initially we need for (|25p a condition on the initial distribution 
function f°(pi): 

J f° Pi d 3 p = 0, (26) 

while from ([2 



f°Pod 3 P=^, (27) 

which is essentially just a normalisation condition. From f)18[> . the finiteness of 
Z{j at r = determines the initial metric, as in [2j, through the equation 



°-%[ /<WA (28) 
a J Po 



(The proof that this determines a?- from f° goes through just as in the earlier 
reference.) Now (j27|) can be seen as just the trace of ([28]) . since at r = 0, 
ib^piPj = bo) 2 - 



S 



For the initial values k®j and Zfj of k%j and Z%j we use dTTJ) and ([20]) evaluated 
at t = to obtain the system: 

— 2^ — a^(fe ) mri £; m?1 = 0, 

7° _l_ i-° 4- f T°\ mntf) _ _ r0 
•"ij ' "'ij T r /ij "'mn ~~ Hj- 

The calculation is precisely as in [2] and we find k®j = Zfj and an equation for 



X%i m {b ) lp {b°) lq k m - 2k %=X%: (29) 



where 



o 167r f f o d3 P 

Xijki = "Ttt / / PiPjPkPi; 



X % = ^ f f\C\ n (b»T m p jPk p m + C k jn (b»n 



+c\ n (b»Y m PlPj p. 



(Po) 2 ' 



which can be thought of as related to the fourth and third moments respectively 
of the initial distribution function, much as (|28|) relates the initial metric to its 
second moment, and (|26p is the vanishing of the first moment. These two 
equations are just as in [2], except with a minor typo corrected. As argued 
there, they give a unique k®- from f°, which is also trace- free, so that Zf, = fe?-. 



2.3 The existence theorem 

In this section, our purpose is to put the equations into a form in which the 
existence theorem of [2j, which was in turn based on work in [9] and [10J, can 
be applied. We shall obtain short-time existence and well-posed-ness for the 
massive Einstein- Vlasov equations with cosmological constant and data at an 
isotropic singularity. The proof is easily simplified to cover existence for the 
limit of massless Einstein- Vlasov. The result is: 

Theorem 2.1 Let G be a 3- dimensional Lie group of some Bianchi type, and 
let (a 1 ) be a basis of left-invariant one-forms on G; let f (pi) be a smooth 
function on the cotangent bundle of G, where pi are the components of p in 
the chosen frame; suppose that f° is compactly supported, supported outside a 
neighbourhood of the origin, and satisfies the constraint 126\): let m be a positive 
constant; then there exists a positive real T and exactly one smooth solution 
(f(t,Pi),9ab) of the space-time Einstein- Vlasov equations on G x (0,T] with an 
isotropic singularity at t = and f(t,pi) — > f°(pi) as t — ► 0. 

We start by transforming (|17p and (|19p . The transformation makes repeated 
use of the following elementary formula: 

1 f 1 ■ 

~(F(t) - F(0)) = / F(sT)ds. (30) 

T JO 
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Introduce new variables Ky and Cij by 

Kij - 
Cij = Z{j — Z^j 



k - - k° 



so that Kij and Cij are initially zero, then (|17p leads to 

—j^Kij = ~(^Cij ~ ^Kij — dijib ) K mn ) 

+ (X°m + K im)(kj n + Kj n )b mn — -b" 1 ™^^ + K mn )(k®j + Kij) 

+P i / nn (k° mn + K mn ) - 8Rij + 8Ar 2 a i3 , (31) 



where, following (|30|) . 



p mn = [\ b mn {k + ^ _ ^mp^O + Kpq ))( ST ) ds . 
JO 



Similarly J20J) leads to 



^Cij — —( Cij K ij (J°)ij nn ^mn) Lij + + (32) 



where 



S$ = -J i tj {sT)d8 

4 2) = -(jfvv)*) (*&+««)■ 

Note that Ijj and J^ 9 can be calculated from (fl4"p . (fT5|) and (pT6|) in terms of 
a, b, k and /. 

Equations (f3Tj) and ([32]) can be written 

^ + ~A^ = G(a, &,«,/), (33) 
ar r 

where u = (k^-, Cij) an d N is a constant matrix constructed from a° and 
(6°)«. In fact iV is exactly as in [2] so that, as there, it is diagonalisable and all 
eigenvalues have positive real part. Suppose it is diagonalised by the change of 
coordinates (y a ) = y = Lu, with a constant matrix L, then (|33p becomes 

dv a 1 

f + -AY = r(a,M,/), (34) 
ar r 

where the eigenvalues X a all have positive real part, and the summation con- 
vention is suspended. The solution of a system like (|34p is dealt with in [9]. It 
can be solved in a form convenient for iteration by the explicit formula 



[ T a x -H a (a,b,y,f)d ( r. 
Jo 
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Equation (|14p is solved by the method of characteristics, but first we use a device 
to deal with the origin in p$ , which is significant for bounding the integrals over 
Pi in (|17p and (|19|) (this is exactly as in [2] which in turn follows |1U|). By 
assumption, f° has compact support and is bounded away from the origin. 
Suppose that f° = for \p\ < B\ and choose B^ < B\ and B% < B2 and a 
smooth function (fi(r) with (j) = for r < B3, eft = 1 for r > B2 and < <j) < 1 
elsewhere. Suppose 6 y is known, and solve 

A p. 

^ = -2cj>(\P\)(b mn P m P n + m 2 T 3 ) 1 ' 2 V l C k ij P k Pi (35) 
for Pi(X;T,pj) subject to 

p i(T;r,pj) =pi. 
Then (I14|) . with a factor c/> on the right-hand-side, is solved by 

f(T,p i ) = f°(P i (0;T, Pj )). 

As long as / has support where (j) = 1 this is the same equation. 

The method of solution for the whole system, as in [2], is now an iteration. 
Given smooth iterates (a n ,b n ,y n , f n ): 

• solve ([35]) using b n and obtain f n +\, this is smooth if b n and f arej 

• solve (fT5]l . (fT6|) and (|34^) with (a n , b n , y n , / n +i) on the right-hand-side for 
smooth (a n+ i,6 n+ i,y n+ i); 

• following the argument in [2] (equations (92)- (97) there), show inductively 
that the iterates are bounded uniformly in n and likewise for the support 
of f n , for some interval in r, say [0, T); 

• now the proof of convergence follows the argument in [2] , the only change 
being an extra term in equation (100) there which becomes 

\y n +i-y n \ < C / {(|a n -a n _i| + |6 n -6 n _i| + |y n -y n _i| + ||/ ri+ i-/ n || 00 )(s) 
Jo 

+ i^J (\ a n ~ a n -i\ + \b n - b n -i\ + \y n - y n -i\ + ||/„+i - fn\\oo)(sp)dp\ }ds; 

• finally uniqueness goes through as before. 



These solutions exist for at least a finite conformal time, and therefore for 
a finite proper time, when they provide data for the existence theorem of Lee 
[5]. Thus for Bianchi types other than type- IX, they will exist forever in proper 
time and have the asymptotic behaviour found by Lee. The metric can be 
rescaled to add the future boundary so, in particular, these solutions expand 
from a bang to a Scri and the conformal metric is extendible through both. 

To obtain a similar result for the massless case, we first need local existence 
with data at the bang. This was done without A in [2]. With A, we seek to 
set m = in the above proof. Note the explicit appearances of m in (I14|) . (|17p 
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and (|19p and recall its implicit appearance in po via (|13p . The first three can 
harmlessly be set to zero (there is an effect on the rate at which the support 
of / spreads, but this is covered by the calculation in [2]). The fourth can be 
set to zero since the support of / is bounded away from the origin, so that 
the integrals over p are finite. Thus the same proof goes through for massless 
Einstein- Vlasov and we have: 

Corollary 2.1 Theorem 0O1 holds for m = 

To find a Scri in the massless case, we next seek to modify the proof of long-time 
existence in [5] to cover m = 0. 

2.4 Long-time existence for massless Einstein- Vlasov with A > 

Long-time existence for some spatially-homogeneous cosmologies whose source 
is massive Einstein- Vlasov with positive A was proved by Lee [3] with the aid of 
a continuation result of Rendall [10] . We now show how this can be extended to 
massless Einstein- Vlasov with positive A. We revert to the space-time metric 
(jlip and proper-time t and by transforming the system (|15 p - (|17p and setting 
m = 0, obtain the system 

j^^Hj = (36) 

±~b» = -b im V n k nm (37) 
at 

~j~kij — 2iRij ~\~ k{ m kj n b 2 kmn)kij ~i~ 2Ag^j 

167T f ~ d 3 P 

J fPiPj^r (38) 

where Rij is as before and 

P° = (b mn p mP n) 1/2 - 

The Vlasov equation becomes 

^ = {b™p m Pn)- ll2 C\ jPm b l i^. (39) 

These equations can be checked against the system in [5] (with some minor 
changes of convention) , and we want to see that the proof of long-time existence 
given there carries over. 

The Hamiltonian constraint is the equation 

C := / fp°d 3 p - WRij + 2A + \{k ij kij - P) = 0, (40) 
yd J 4 

and the momentum constraint follows from (1251) as 



Ci := t-^d = ^| J f Pt d 3 p - 2b mn (C\J nj + C j mj k ni ) = 0. (41) 
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Theorem 2.2 Given data (dij,b lJ , fey, /) at some time to, with k > and sat- 
isfying the constraints |^0| ) and |^ i| ), i/ten for Bianchi types other than type-IX, 
there exists a unique solution of the Einstein equations with this data, which 
furthermore exists for all t > to. The solution admits a Scri. 

The proof is very similar to that in [5] . From the trace of (|38|) , with the aid 
of (|4"0]1 . we obtain 



it 

Decompose as 



d ± = 2A - - ^ I ffa*p. (42) 



'a 



kij — 0~ij -\- ^ kdij , 



so that o~ij is trace-free, then (|40|) can be rearranged to give 
16tt 



'a 



J fp°d 3 p - b i! R u + \o^G i3 = -2A + (43) 



Following [H], we restrict to spatially-homogeneous metrics not of Bianchi type- 
IX, then these have b t3 Rij < 0, so that the left-hand-side in (|43p is non- negative. 
Putting this with (|4"2"j) we obtain 

d£; . 1-9 
— < 2A - -k 2 < 0, 
~ 6 ~ 

for all but type-IX. Now k is initially positive and we can integrate this last 
equation to find 

6H<k< 6H y_ e _ 2 J y (44) 



where H = A/3, which is the 7 of [5] . This bounds k above and below for 
positive t and forces k — ► as i — >• 00. From f|36|) 

-loga = fc, 

so that log d is also bounded above and below for positive t and tends to a linear 
function of t. From (I43p . 



1 7 v * lr 2 8Ae~ 2Ht 

~o~: o~~ < -2A + -k < 



3 - 1 6 - (l-e- 2 ^) 2 ' 

so that <t. j is bounded and tends to zero, and we may follow the argument in 
[5] to find the asymptotic form of the metric. The solutions exist for ever, are 
complete in the future and can be rescaled to add future infinity. 

Putting this with the Corollary [2Tj we obtain solutions of massless Einstein- 
Vlasov with positive A which expand forever from an isotropic singularity. From 
the work of Lee, [5], it follows that the distribution function has a limit at Scri, 
but there seems no reason to expect that ([26]) of ([28]) will hold at Scri. If 
this is so, then for this matter model bang and Scri have different conformal 
properties. 
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3 Radiation with A 



In this section, the object is to prove existence of examples of spatially-homogeneous 
cosmological models, whose source is a perfect fluid with the radiation equation 
of state and a cosmological constant, with data given at an isotropic singularity. 
We also show that, for positive A, some of these expand for ever (in proper time) 
and have a smooth Scri. We present a series of increasingly complex examples, 
starting with a simple case. 

3.1 The example of FRW 

The FRW space-time metric is 

g = dt 2 - R 2 (t)dal 

where R(t) is the scale factor and do\ stands for the metric of a 3-dimensional 
space of constant curvature which is positive, negative or zero according as k 
is 1, —1 or 0. 

With a perfect fluid source and the radiation equation of state, p = p/3, the 
Einstein field equations reduce to the conservation equation, which integrates 
to give 

pR 4 = m, (45) 

where m is a (positive) constant of integration, together with the Friedman 
equation 

te + k = ^-pR 2 + |i? 2 . (46) 

Here the over-dot indicates differentiation with respect to proper time t. We 
introduce conformal time r by 

dr = dt/R(t) (47) 

and denote differentiation with respect to it by a prime. Then using (|45p in 
([MD and setting = l we find 

(R'f = j~kR 2 + V. (48) 

This equation has a unique solution with R = 0, R' > at r = 0, corresponding 
to an initial singularity. If we assume that either k is negative or zero or that 
k = 1 but mA > 9/4 then this solution blows up in finite r-time, at say r = Tp 
given by 

dR 

(™ _ kR? + |i? 4 )V2 ' 

In fact R has a simple pole R ~ {H(rp — t)) -1 where, as in the previous 
section, H = y/A/3. Thus, by (I47j) . in terms of proper time R ~ H~ 1 e Ht and 
Tp corresponds to future infinity or Scri. 
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Equation (|48p also has a striking symmetry: if we replace R by R = R 
we obtain the same equation but with m and A interchanged: 

(Elf = A _ k R2 + ™£4 

K J 3 3 

Evidently this symmetry interchanges Scri and the initial singularity, while also 
interchanging the constants m and A. The conformal metric extends smoothly 
through both ends, but this is perhaps unsurprising as the metric is conformally 
flat. 

In subsequent examples we find that the conformal extension property is 
preserved, but not the symmetry interchanging Scri and bang, which is dis- 
rupted by the presence of conformal curvature. 

3.2 Bianchi-III and Kantowski-Sachs models 

These have a space-time metric which can be parametrised as follows: 

g = dt 2 - {RS 2 fdz 2 - (RS' 1 ) 2 '{d6 2 + f k {9) 2 d(j) 2 ). (49) 

Here k = +1 and f\ = sin# for the Kantowski-Sachs metric, and k = — 1 and 
/_i = sinh# for the Bianchi type-Ill metric. 

This choice of parametrisation allows the introduction of conformal time r 
just as in (|47l) and the solution of the conservation equation just as in (145 j) . 
The conformally rescaled metric 

g ■= RT 2 g = dr 2 - S 4 dz 2 - S~ 2 {d8 2 + / fe (^) 2 # 2 ) 

is evidently smooth wherever S is smooth and non-zero, and we shall see ex- 
amples where this includes the initial singularity and Scri. 
The Einstein equations can be written as the system 

4 +6 H- 3 (l) - ks2 - ° < 5 »> 

Till / Ol \ 2 

3^ + 3(^j +kS 2 -2KR 2 = (51) 

where these have been simplified with the aid of the Hamiltonian constraint, 
which in turn can be written as 

3 (l) 2 - 3 (l) 2 - ks2 + S +Afi2 < 52 » 

We shall see below that these have solutions with initial data (R, R' , S, S') = 
(0, ym/3, Sq, 0) at (say) r = 0. Suppose this is true, then given A, this is a 
2-parameter family expanding from an isotropic cosmological singularity. The 
Kantowski-Sachs metric could recollapse but the Bianchi type-Ill metric cannot, 
as we see as follows: introduce 

5 /x 2 



Q = 3[-) -kS 2 , (53) 
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which is manifestly non-negative for k = — 1. Then by (|50p 

so that, with the chosen initial conditions, 

< Q < Si 

Therefore S and ^ are bounded for all time. Write ([52]) as 

The right-hand-side is strictly positive, so that R' , positive initially, never van- 
ishes. Also ^ is bounded as long as R is. Thus solutions exist until R diverges. 
This happens after finite conformal time, as we see by comparing R with the 
solution L of the equation 

3(L') 2 = m + AL 4 , L(0) = 0, 1/(0) > 0. 

Then R> L but L diverges in a time 

dL 

(m + AL 4 ) 1 / 2 ' 

which is finite. We shall return to the asymptotic form after proving existence 
of solutions with data as claimed. For this, we put the system of Einstein 
equations into a first-order Fuchsian form, like (|33|) . Set 

R = T e u , S = e s , (55) 

when ([SU]) and (|5"Tj) become 

£' = Z 
[/' = VF 

r 3 

W = --^-Z 2 -W 2 --e 2S + — r 2 e 2U 
t 3 3 

while the Hamiltonian constraint becomes 

m = e 2U (3 + 6tW + t 2 (3W 2 - 3Z 2 + ke 2E - Ar 2 e 2U )) , 

which can be interpreted just as the statement that the right-hand-side is con- 
stant, which evaluation at r = shows to be positive. (In this setting, the 
momentum constraint is vacuously satisfied.) 

This is the desired first-order Fuchsian form and satisfies the conditions for 
Theorem 1 of [9j, so that it has smooth solutions given data (£, U, Z,W) = 
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(So, Co, 0, 0) at r = 0. Clearly the two free parameters correspond to m and So 
as claimed. 

It is straightforward to show inductively that as power series about r = 
U and £ have only even powers of r so that S is even while R is odd. In 
particular the initial singularity is umbilic in the conformal metric. We try the 
corresponding expansion at Scri: if R diverges at say r = Tp then introduce 
R = R^ 1 so that (f54l) becomes 

3(R') 2 = A + QR 2 + mR 4 , (56) 



and, as t — ► Tp, R — ► and i?' — ► — yA/3 = — H. Thus i? has a power series 
(strictly, a Laurent series) beginning with a first-order pole. We can write (|50p 
as 

#2^' = = _^ R 2 S 2 (57) 



5^3 3 
in the Bianchi type-III case, so that 

,S' 1 



R 2 ^ = -\[ Q R 2 {*)S\*)do, 



from which, given the asymptotic behaviour of R just found, — ► and 
<5 — ► 5*00 say, as t ^ Tp. Now differentiate ([56]) again to find 



3R" = Q + 6 — i? + 2mi? 3 , 




so that i?" — > as r — ► Tp, and differentiate again to find 

it — > ---H ^oo- 

From the integral solution of (|57p we find 5" and «S"" finite at Tp and, if we set 
f = Tp — t, then power series begin 

R = C(f + ^|f 3 + -..) 
i2 = H- 1 f- 1 (l-^f 2 + ---) 
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5; 



3 



S = S 00 + -fff 2 + S 3 f 3 + ---. 

Here Sqo and 53 are fixed by initial conditions but not in a way we can find 
explicitly. From the point of view of Scri, they are free data. Scri is located 
at f = and is umbilic in the conformal metric, just as the initial singularity 
was. We have done enough to show that the metric is at least C 3 there, and 
inductively we could get smoothness, but S is not necessarily even about Scri. 

The fact that Soo is free is equivalent to saying that the metric of Scri is 
free data but S3 corresponds to part of the Weyl tensor (more accurately, to 
the time-derivative of the magnetic part of the Weyl tensor). See in this con- 
nection the discussion of vacuum solutions with data at Scri in [IT], and of the 
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asymptotic behaviour of solutions in [12J: the time-derivative of the magnetic 
part of the Weyl tensor is free data at Scri but not at the initial singularity. Of 
course, for the particular solutions considered here, which expand from initial 
isotropic singularities, the Weyl tensor at Scri is determined by the data at the 
bang, and it remains a possibility that S3 is obliged to vanish at Scri. However, 
from the point of view of initial value problems, there is more free data at Scri 
than at the bang and if one evolves back from Scri with arbitrary choices of the 
data there then one cannot expect to arrive at an isotropic singularity at the 
beginning. 

Another way to see this last fact is to write the Einstein equations as a 
first-order system about Scri by introducing R = f _1 C7 and S = e s in place 
of (|55p . The resulting first-order system is Fuchsian but does not satisfy the 
conditions for Theorem 1 of [9] (since it has to allow extra data at 0(f 3 ) in 
£). This is one explanation for the absence of the inversion symmetry found 
for FRW. 

3.3 Bianchi class A models 

With fluid flow orthogonal to the surfaces of homogeneity, there is no loss of 
generality in assuming that the metric is diagonal in the invariant basis. Thus 
we may write the space-time metric as 

g = dt 2 - R 2 (e 2a a 2 + e 2/ V 2 2 + e 2 V 3 2 ), (58) 

where the Oi are a basis of left-invariant one- forms for the isometry group. For 
Bianchi class A the one-forms can be assumed to satisfy the system 

da 1 = ni<72 A 03 
da 2 = n20"3 A o\ 
do 3 = n^a\ A 02 

where the rii are chosen from (1,-1,0) (and inequivalent possibilities lead to 
different Bianchi types, see e.g. [8]). The parametrisation of the metric is 
chosen so that R G is the determinant of the physical spatial metric, so that 

a + /3 + 7 = 0. 

We can, as in the previous two cases, solve the conservation equation to find 

pR A = m = constant 

and we may introduce conformal time r as in (I47p . Then the conformally- 
rescaled metric 

g : = R-2~ g = dr 2 - (e 2a a 2 + e 2 % 2 + e 2 V 3 2 ) 

is smooth wherever (a, /3, 7) are. It gives the metric of the bang where R = 
and the metric of Scri where R diverges. 
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It will be convenient to write 

iVi=n ie 2Q , jV 2 = n 2 e 2/3 , N 3 = n 3 e 2 \ 
when we may write the Einstein equations at once in terms of t as 

a" + 2^-a' + A = 0, (59) 
R 

and two cyclic permutations of this for (5 and 7. Here 

A = ~(2N? - N 2 2 - N 3 2 + 2N 2 N 3 - iViiV 2 - JV3JV1), 

and we write B and C for the two cyclic permutations of this. Together with 
(|59p we need the Hamiltonian constraint which is 

|) 2 = ^((«0 2 + (/?') 2 + (7?) + ^ + Ai? 2 - F (60) 
where 

F = -\{N 2 + N 2 2 + iV 3 2 - 2N 2 N 3 - 2N 3 Ni - 2NiN 2 ), (61) 

which is proportional to the Ricci scalar of the spatial metric. (Again the 
momentum constraint is vacuously satisfied.) By differentiating (|60l) we can 
obtain a second-order equation for R which can be written as 

^ = -\{W) 2 + (P') 2 + (if) ~\f+ IaR 2 - (62) 

Our first aim is to show that solutions exist with isotropic singularities. As 
in the previous subsection, this is accomplished by putting the equations into 
Fuchsian form. We introduce U by R = re u as in (|55p and then put the 
equations into first-order form by introducing variables (X, Y, Z, W) according 
to 

a = X 

= Y 

1 = Z 

u' = w 

Then (j59j) and its permutations become 
X' = —X-2WX-A 

T 

Y' = --Y-2WY-B 

T 

Z' = — Z — 2WZ — C 

T 

W = --W-W 2 --(X 2 +Y 2 + Z 2 )--F + —r 2 e 2U . 
t 6 3 3 
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The combined system of eight equations is in Fuchsian form, and satisfies the 
assumptions of Theorem 1 of [9|. Thus solutions exist in a neighbourhood of 
r = with (a, [3, 7) and U freely specifiable there and (X, Y, Z, W) all zero 
there. 

The Hamiltonian constraint (|60|) can be written as 

m = e 2U U + 6tW + t 2 (3W 2 - ^(X 2 +Y 2 + Z 2 ) + F - At 2 e 2U ) 

when it can be read as relating the (positive) constant of integration m to the 
initial value of U. 

To study the behaviour at late times, we restrict to the Bianchi class A 
metrics other than type-IX. Recall that type-IX has (^1,^2,^3) = (1, 1, 1) and 
F in (|6ip can have either sign. For all the other types, F is non-positive. 
Thus from ([62]) R" is positive, therefore so is R' and solutions expand until R 
diverges. There is a finite conformal-time blow up in R, at tf say, since (I60p 
gives 

3(R') 2 >m + AR 4 

and we can argue just as in subsection 3.2. Again R will have a pole in t: 
R ~ (H(t f - t))~ x . Analogously to flM}, we introduce 

Q = X 2 + Y 2 + Z 2 - IF 

then 

Q> = -^(x 2 + Y 2 + Z 2 ), (63) 
R 

so that Q decreases and, for all time, 

< Q < Q(0) = -6F(0). 

Thus X, Y and Z are bounded for all time, so a, j3 and 7 grow at worse linearly 
in r and, since there is only a finite amount of r-time, are also bounded and 
have limits at rp. These finite limits define the metric of Scri. From (|59j) we 
obtain 

R 2 X := R 2 a' = - [ R 2 (a)A(a)da, (64) 
Jo 

so that a' — > as r — > Tp. Thus X, Y and Z vanish at Scri and Scri is umbilic 
in the rescaled metric. For more differentiability, we mimic the calculation in 
the previous section. First set R = R^ 1 so that ([60]) becomes 

3(R') 2 = A + QR 2 + mR 4 , 

and then with this and (163j) and (]64h we can obtain power series for a, (3, 7 
and R. The rescaled metric will be smooth at Scri. 

As in the previous subsection, R is odd in r and (a, [3, 7) are even about 
the bang, but this doesn't necessarily hold at Scri, where the terms of order 
f 3 in (a, /5, 7) are not determined in the power series. Scri and the bang are 
both umbilic in the conformal metric, but are not equivalent, and in particular 
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we don't have the symmetry interchanging them found for FRW. We do have 
solutions of the kind considered by Penrose [7J, where the conformal metric is 
extendible through the bang and through Scri, but, at least with this matter 
model, it seems that we can't use Scri as a new bang, that is as the initial data 
surface for another cycle of expansion, without having the conformal metric not 
smooth, i.e. the normal-derivative of the magnetic Weyl tensor discontinuous. 
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